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On the spin distributions of ACDM haloes 
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Abstract We used merger trees realizations, predicted 
by the extended Press-Schechter theory, in order to 
study the growth of angular momentum of dark matter 
haloes. Our results showed that: 

1) The spin parameter A' resulting from the above 
method, is an increasing function of the present day 
mass of the halo. The mean value of A' varies from 
0.0343 to 0.0484 for haloes with present day masses in 
the range of 10^h~iA//o to lO^h-^Mo. 

2) The distribution of A' is close to a log-normal , but, as 
it is already found in the results of N-body simulations, 
the match is not satisfactory at the tails of the distribu- 
tion. A new analytical formula that approximates the 
results much more satisfactorily is presented. 

3) The distribution of the values of A' depends only 
weakly on the rcdshift. 

4) The spin parameter of an halo depends on the num- 
ber of recent major mergers. Specifically the spin pa- 
rameter is an increasing function of this number. 

Keywords galaxies: halos - formation -structure; 
methods: numerical -analytical; cosmology: dark mat- 
ter 



1 Introduction 

There are two more likely pictures regarding the growth 
of angular momentum in dark matter haloes. 
The first is that galactic haloes acquired their angu- 
lar momentum from the tidal torques of the surround- 
ing matter. This is an old idea of Hoyle (1949) that 
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has been investigated in a large number of studies (e.g. 
Peebles 1969; Doroshkevich 1970; Efstathiou & Jones 
1979, Barnes & Efstathiou 1987, White 1984; Voglis & 
Hiotehs 1989, Warren ct al. 1992, Stcinmetz & Bartel- 
mann 1995). The results of the above studies, analytical 
and numerical, show that the spin parameter A, intro- 
duced by Peebles (1969) and defined by the relation 
A = Jv/rE]/GA/5/2, has an average value of about 
0.05-0.07, where J is the modulus of the spin of halo, 
M, E are the mass and the energy respectively and G 
is the gravitational constant. According to the above 
picture, haloes acquire their angular momentum dur- 
ing their linear stage of their evolution because during 
this stage they have large linear sizes and thus the en- 
vironment is capable to affect their evolution by tidal 
torques. Since their expansion is decelerating their rel- 
ative linear size becomes smaller and the affection by 
the environment becoms less significant. The moment 
of the maximum expansion is in practice the end of 
the epoch of growth of angular momentum. At latter 
times, the halo evolves as a dynamically isolated sys- 
tem. Steinmetz & Bartelmann (1995) showed that the 
dependence of the probability distribution of A on the 
density parameter of the model Universe as well as on 
the variance of the density contrast field is very weak. 
Only a marginal tendency for A is found to be larger 
for late-forming objects in an open Universe. 
The second picture is closely related to the hierarchi- 
cal clustering scenario of cold dark matter (CDM; Blu- 
menthal et al. 1986). According to this scenario, struc- 
tures in the Universe grow from small initially Gaussian 
density perturbations that progressively detach from 
the general expansion, reach a maximum radius and 
then collapse to form bound objects. Larger haloes 
are formed hierarchically by mergers between smaller 
ones, called progenitors. The buildup of angular mo- 
mentum is a random walk process associated with the 
mass assembly history of the halo's major progenitor. 
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The main role of tidal torques in this picture is to pro- 
duce the random tangential velocities of merging pro- 
genitors. 

The above two pictures of formation are usually stud- 
ied by two different kinds of methods. The first kind 
is the N-body simulations that are able to follow the 
evolution of a large number of particles under the in- 
fluence of the mutual gravity, from initial conditions to 
the present epoch. The second kind consists of semi- 
analytical methods. Among them, the Press-Schechter 
(PS) approach and its extensions (EPS) are of great 
interest since they allow to compute mass functions 
(Press & Schechter 1974; Bond et al. 1991) to ap- 
proximate merging histories (Lacey & Cole 1993, Bower 
1991, Sheth & Lemson 1999b) and to estimate the spa- 
tial clustering of dark matter haloes (Mo & White 1996; 
Catclan et al. 1998, Sheth & Lemson 1999a). 
Recently large cosmological N-body simulations have 
been performed in order to study the angular momen- 
tum of dark matter haloes in ACDM models of the Uni- 
verse (e.g. Bullock 2001, Kasun & Evrard 2005, Bailin 
& Steinmetz 2005, Avila-Reese et al. 2005, Gottlober 
& Turchaninov 2006). 

Additionally, semi-analytical methods like merging his- 
tories resulting from EPS methods, have been used for 
similar studies (Vitvitska et. al 2002, Mailer et. al 
2002). 

In this paper, we use such merging histories based on 
EPS approximations to study the distribution of spins. 
The paper is organized as follows: In Sect. 2, basic equa- 
tions are summarized. In Sect. 3, we present our results 
while a discussion is given in Sect. 4. 



2 Construction of merger-trees and acquisition 
of angular momentum 

According to the hierarchical scenarios of structure for- 
mation, a region collapses at time t if its overdensity at 
that time exceeds some threshold. The linear extrapo- 
lation of this threshold up to the present time is called 
a barrier, B. A likely form of this barrier is: 



B{S, t) = VoS; [1 + l3{S/aS,)' 



(1) 



In the above Eq. a, (3 and 7 are constants, S** = 
S^,{t) = S'^{t) where 6c{t) is the linear extrapolation 
up to the present day of the initial overdensity of a 
spherically symmetric region that collapsed at time t. 
Additionally, S = cr^(M), where cr^(M) is the present 
day mass dispersion on comoving scale containing mass 
M. S depends on the assumed power spectrum. The 
spherical collapse model has a barrier that does not 
depend on the mass (eg. Lacey & Cole 1993). For 




Fig. 1 The figure shows the conditions at the onset of the 
merger of two haloes with masses mi and m2 and virial 
radii Ri and R2 respectively with mi > m2 and Ri > i?2. 
The centers of their masses are indicated by CMi and CM2 
and their distance is r = max{Ri, R2) — Ri- The position 
vectors of tlieir centers of mass are ri and r2 respectively 
so their relative position is CM12 ~ r2 — ri. The vectors of 
the velocities of their center of mass are vi and V2 and the 
vector of their relative velocity is v^ei = V2 — vi. Haloes 
merge if they approach each other, that is the condition 
Vrei ■ CM12 < is fulfilled. See text for more details. 



this model, the values of the parameters are a = 1 and 
P = 0. The ehipsoidal collapse model (EC) (Sheth, 
Mo & Tormen 2001) has a barrier that depends on the 
mass (moving barrier). The values of the parameters 
are a = 0.707,/3 = 0.485, 7 = 0.615 and are adopted 
either from the dynamics of ellipsoidal collapse or from 
fits to the results of N-body simulations. 
Sheth & Tormen (2002) showed that given a mass ele- 
ment -that is a part of a halo of mass Mq at time Iq- 
the probability that at earlier time t this mass element 
was a part of a smaller halo with mass Af , is given by 
the equation: 



f{S,t/So,to)dS = 



1 |r(^,t/^o,t o)| 

2^ (A5)3/2 



■ exp 



(AB)2 
2A5 



dSl 



(2) 



where AS S - Sq and AB = B{S, t) - B{So, h) with 

S^S{M),So^S{Mo). 

The function T is given by: 



T{S,t/So,to) = B{S,t)-B{So,to)+Y^ 



n=l 



(3) 



BiS,t).l 
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Fig. 2 Spin parameter A' for the most massive progenitor 
at scale factor a for the cases 1, 2 and 3. SoUd, dashed and 
dotted lines correspond to the cases 1, 2 and 3 respectively. 
The evolution of A' is characterized by sharp increases and 
decreases due to mergers. 



Recently, Zhang & Hui (2006) derived first crossing dis- 
tributions of random walks with a moving barrier of an 
arbitrary shape. They showed that this distribution 
satisfies an integral equation that can be solved by a 
simple matrix inversion. They compared the predic- 
tions of their exact numerical solution with those of 
approximation given by cq.3 and found a very good 
agreement. This shows that eq. 3 works well for mildly 
non linear barriers as that given by eq.l above. 
Eq. 2. can obviously predict the unconditional mass 
probability, t), which is simply the probability that 
a mass element is at time t a part of a halo of mass M, 
by setting = 0, and B{So,to) = 0. We note that the 
quantity Sf{S,t) is a function of the variable v alone, 
where = 5c[t)/a{M). 5c and a evolve with time in 
the same way, thus the quantity Sf{S,t) is indepen- 
dent of time. Setting 25/(5', t) = vfiv) one obtains 
the so-called multiplicity function f{i^). The multiplic- 
ity function is the distribution of first crossings of a 
barrier B^iy) (that is why the shape of the barrier in- 
fiuences the form of the multiplicity function) , by inde- 
pendent uncorrelated Brownian random walks (Bond ct 
al. 1991). The multiplicity function is related to the 
comovimg number density of haloes of mass M at time 
t, N{M,t), by the relation. 



r, ^ ,./,r ^dlnM 
t^fH = ^TN(M,t)— 



(4) 



model of the Universe at time t. 

Using a barrier of the form of Eq. 1 in the unconditional 
mass probability, one finds for f{iy) the expression: 

/(z.) = y^[l+(3{a,^^y^g{j)] exp {-0.5a,^^[l + f3{a,y^r^]^)l 

(5) 



where 



5(7) =1 1-7 + 



7(7 - 1) 



2' 



7(7- 1) ■ ■ ■ (7-^4) 
5! 



(6) 



that results from the excursion set approach (Bond et 
al. 1991). In the above Eq., pb{t) is the density of the 



Recent comparisons show that the use of EC model im- 
proves the agreement between the results of EPS meth- 
ods and those of N-body simulations. For example, 
Yahagi et al. (2004) show that the multiplicity func- 
tion resulting from N-body simulations is far from the 
predictions of spherical model while it shows an excel- 
lent agreement with the results of the EC model. On 
the other hand, Lin et al. (2003) compared the distri- 
bution of formation times of haloes formed in N-body 
simulations with the formation times of haloes formed 
in terms of the spherical collapse model of the EPS the- 
ory. They found that N-body simulations give smaller 
formation times (i.e. higher redshifts). Hiotelis & Del 
Popolo (2006) showed that using the EC model, for- 
mation times are shifted to smaller values than those 
predicted by a spherical collapse model. Additionally, 
EC model, combined with the 'stable clustering' hy- 
pothesis, was used in order to study density profiles of 
dark matter haloes (Hiotelis 2006). The resulting den- 
sity profiles at the central regions of haloes have the 
interesting feature to be closer to the results of obser- 
vations than the results of N-body simulations. Conse- 
quently, the EC model is a significant improvement of 
the spherical model and therefore we use this model for 
our calculations. 

We assume a number of haloes with the same present 
day mass Mq -at present epoch to- and we study their 
past using merger-trees by finding their progenitors - 
haloes that merged and formed the present day haloes- 
at previous times. The procedure for a single halo is as 
follows: A new time t < is chosen. Then, a value AS 
is chosen from the desired distribution given by Eq.2. 
The mass Mp of a progenitor is found by solving for Mp 
the equation AS = S{Mp) - S{Mq). If the mass left to 
be resolved Mq — Mp is large enough, the above proce- 
dure is repeated so a distribution of the progenitors of 
the halo is created at t. If the mass left to be resolved 
-that equals to Mq minus the sum of the masses of its 
progenitors- is less than a threshold then, we proceed 
to the next time analyzing with the same procedure the 
mass of each progenitor. The most massive progenitor 
at t is considered as the mass of the initial halo at that 
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time. 

A complete description of the above numerical method 
is given in Hiotclis & Del Popolo (2006). The algo- 
rithm - known as N-branch merger-tree- is based on 
the pioneered works of Lace & Cole (1993), Somcrvillc 
& KoUat (1999) and van den Bosch (2002). 
In our calculations, we used a flat model for the Uni- 
verse with present day density parameters Qm.o ~ 0.3 
and r^A.o = " ^ the cosmological con- 

stant and Hq is the present day value of Hubble's con- 
stant. We used the value Hq = lOOhKMs-^Mpc"! and 
a system of units with rriunit = lO^^Moft,"^, Tunit = 
l/i~^Mpc and a gravitational constant G = 1. At this 
system of units, Ho/Hunit = 1.5276. 
As regards the power spectrum, we used the ACDM 
form proposed by Smith et al. (1998). The power spec- 
trum is smoothed using the top-hat window function 
and is normalized for ag = a{R = 8/i~^Mpc) = 1. 
A merger-tree gives the complete history of the halo. 
After its construction, we know all the progenitors of 
a halo with present-day, at time t — to, mass Mq at a 
previous time ti, all progenitors of every progenitor at 
ti at time t2 < ti etc. This procedure is repeated up 
to a level n of resolution that corresponds to a time i„ . 
As regards the choice of the time-step we used the pro- 
cedure described in Hiotelis & Del Popolo 2006 that is: 
The equation 5c{anew) = Aw + (5c(ao;d), where aou = 1 
at the beginning of the construction, is solved for anew , 
that is the new value of the scale factor. We used a 
constant value of Acj = 0.1 but tests with smaller val- 
ues of Acj showed no differences in the results. The 
procedure stops for a < amin = 0.05. Then, we start to 
merge haloes present at time t„ to create the angular 
momentum history of those haloes that are present to 
the time level tn-i- We assume that two haloes with 
virial masses toi and TO2 and virial radii ri and r2 merge 
when the following conditions are fulfilled: 

a) They approach each other. 

b) Their relative energy, given by = ^fj.v'^ — '^"^^"'^ , 



is negative, where fj, 



r and v are the distance 



and the relative velocity of their canters of masses, 
c) Their distance r is equal to the the maximum of ri 
and r2. 

After such a merge a new halo is created with mass 
'7^1,2 = 'm-i + '7i2 and spin Si, 2 given by 



Si, 2 — Si 



-'orb,1.2 



(7) 



where Si and S2 are the spins of the two haloes and 
Lorb,i.2 is their orbital angular momentum given by 



-iorb.1,2 



mim2 

mi + ni2 



(r X v) 



(8) 



r and v are the vectors of relative position and veloc- 
ity of their center of mass respectively. A halo which 



has suffered no merger up to a time t has no spin and 
consequently all haloes at time t„ have no spin. The 
virial r radius of an halo at scale factor a is related to 
its virial mass m by the relation: 



r{a) = 

where 
X(a)E 



2Gm{a) 



1/3 



~(a)rim(a)-ff^(a)- 



(9) 



(10) 



r2m(a) and H{a) are the density parameter and the 
Hubble's constant at scale factor a, respectively. For 
Aitir we used the expression given in Bryan & Norman 
(1998) Ay,r{a) ~ (187r2 - 82a; - 39x^)/n,n{a) where 
X = I — flm{a). 

The construction of a merger-tree does not requires 
or predicts any information about the velocity field of 
merging haloes. Since, in our case, the purpose is to 
study the growth of angular momentum during a pro- 
cess of subsequent mergers, we need a model for the 
description of the velocity field. So, at first, we used 
an arbitrary, but reasonable, model, that is described 
in details below, satisfying the conditions a) b) and c) 
set above. Additionally in section 3. we also refer to a 
model for the velocity field that is consistent with the 
results of N-body simulations. Both models gave simi- 
lar results, but the question which model describes the 
velocity field best is open and under investigation. 
The whole procedure of merging the haloes present in 
the merger-tree model follows: 

Let that at level I, a set k haloes with virial masses 
m,i,m,2---m,k, virial radii '"i,2---7';,fc and spins 
Si,i, Si,2...S;,fc consists of all the progenitors of an halo 
with virial mass m;_i.i and virial radius ri_i,i at level 
I — 1. The merger procedure is as follows: Two pro- 
genitors mi,i and to;, 2 merge. First r = maa;(n,i, ri,2) 



and 



are calculated. Then, the maxi- 



»Tii,i-l-"i;,2 

mum relative velocity, that satisfies the condition of 
negative total orbital energy, is found by the relation 

Vmax = ( ^'^""■i"'.^ )l/2^ Thcu, thc moduluS, Vrel; of 

the relative velocity vector v^ez of the two haloes is 
picked by a Gaussian distribution with mean value 



c/2 and <t = (l/3)w„ 



The two 



components and Vy are found using uniform dis- 
tributions in the range [—Vrei,Vrei\- If the condition 
t'x + — ''^rei '^^ Satisfied, the third component if 
found by = ± ^ v^^^ — [v"^ + ''^y ) j where the sign is 
chosen randomly. If the above inequality is not satis- 
fied, new values of Vx and Vy are chosen and the proce- 
dure is repeated. The components [x, y, z)of the rela- 
tive position vector r, with modulus r = max(ri^i, r;,2) 
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arc defined choosing x and y by a uniform distribu- 
tion in the range [— r, r] and then, if the condition 
+ < is fulfiUed , the component z is defined by 
hy z = ± -^/r^ — .T^ — . The condition of approach- 
ing, r • Vrad < 0, where Vrad is the radial component of 
the relative velocity, is checked. If the condition is not 
fulfilled we go back to choose new velocity components, 
otherwise we continue by finding the orbital angular 
momentum, according to (8) and the spin of the newly 
formed halo according to (7). The mass of this halo is 
mi,i +TOi,2, and its virial radius is defined by (9) where 
a is the scale factor of the Universe at level I. The 
number of k haloes of the set is now reduced by one. 
The procedure is repeated until the number of haloes 
becomes one. Thus, the angular momentum of a new 
halo, present at time level l — l, which had at time level 
I the above k progenitors, is found. The procedure is 
repeated for all the sets of progenitors at level I and so 
the angular momentum for every halo of level ? — 1 is 
calculated. This new set of haloes consists of the pro- 
genitors of haloes at level 1 — 2. The same procedure is 
repeated for level ^ — 1 to create the haloes at level I — 2 
and so on. The procedure ends with the formation of 
a single halo of mass mo at level that represents the 
present age to of the Universe. Fig.l shows two haloes 
at the onset of their merger. 

As a measure of the angular momentum, we used the 
parameter A' given by the (see Bullock et al. 2001; 
Dekel et al. 2000) 

A' - 7^ (11) 

where S equals to the spin, m is the virial mass, r the 
virial radius of the halo and Vc = \/ Gm/r is the circular 
velocity at distance r. The parameter A' is easier to be 
measured, than A, not only in simulations but in semi- 
analytical methods as the merger-tree method used in 
this paper. This happens because the total energy of 
the halo is not required. Instead, the calculation of 
A requires a known density profile p for the halo, the 
calculation of the potential cj) by the expression 

$(r) -47rG[- / p(r')r'^dr' + / p(r')r'dr'], (12) 

^ Jo Jr 

the calculation of the potential energy W 

W J p(r)$(r)47rr2dr (13) 

and finally, assuming virial equilibrium, the derivation 
of the total energy hy E = ^ly. Simplified forms for 
specific density profiles can be found in Sheth et. al 
2001 



It is noticed that the spin parameters A and A' are ap- 
proximately equal for typical Navarro et al. 1997 haloes 
(Bullock et al. 2001). 

3 Results 

We studied five cases for haloes of different present day 
masses ttiq. In our system of units, nio takes the values 
10-2, 0.1,1., 10 and 100 for the cases 1,2,3,4 and 5, 
respectively. For every case, we produced a number of 
Nres = 20000 — 60000 realizations. Studying the pro- 
genitors of a halo with present day mass toq at some 
time tt. it is likely to find a number of them with zero 
angular momenta. This is a consequence of the way 
of the construction of the merger-tree. These progen- 
itors are haloes that have not suffered any merger for 
times t < and thus either they have not increase 
their masses at all or they have accreted only small 
amounts of matter. Notice that the growth of angu- 
lar momentum resulting by the accretion of amounts of 
matter that are below a critical value -see Hiotelis & 
Del Popolo (2006), for the details of the construction of 
merger trees- is not taken into account in our scheme. 
The distributions presented below are predicted by tak- 
ing into account only haloes that have no zero angular 
momentum. 

We note here that the growth of angular momentum 
in the above presented picture is not so smooth as in 
the case of the tidal torque theory where the angular 
momentum is a simple increasing function of time (eg. 
Barnes & Efstathiou 1987, Voglis & Hiotelis 1989). In 
the picture studied in this paper, angular momentum 
increases or decreases in a complicated way. As an ex- 
ample. Fig. 2 shows the evolution of A' as a function of 
the scale factor a for randomly selected histories from 
the cases 1,2 and 3. This picture is characterized by 
sharp increases and decreases of the spin parameter 
due to mergers. This is a common behavior present 
in similar studies (e.g. Vitvitska et. al 2002). In the 
following we study four important characteristics, of 
the above picture of the growth of the angular momen- 
tum, namely: 1) The form of the distribution of A'. 2) 
The dependence of this distribution on the present day 
mass of the halo. 3) The role of major mergers and 
their affection on the magnitude and the distribution 
of A', and 4) the dependence of the distribution on the 
redshift. 
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Fig. 3 The present-day, at 2 = 0, distributions of A' for the cases CI and C5 are given by squares at the left and the right 
snapshots, respectively. Dashed lines are the best least-squared fits by a log-normal distribution given in text by eq. 14. 
Solid lines are also the best least-squared by a distribution given by eq. 15. It is clear that the fits by eq. 15 approximate 
better the results than the log-normal distribution. In case C5, this is very clear at the right tail of the distribution. 



It has been proposed in the literature that the re- 
sulting distribution of A' is well approximated by a log- 
normal distribution with probability density given by 



0(A') 



exp 



ln2(A7A'^ 



2al, 



(14) 



(e.g. van den Bosch 1998; Gardner 2001; Bailin & 
Steinmetz 2005). This is a Gaussian in In A'. The ex- 
pectation value for In A' is < In A' >= In A' and it peaks 
at Apg^j. = \' exp{—a\,). A recent analysis of large cos- 
mological simulation was performed by Bett et al.2007. 
These authors used the results of the Millenium simu- 
lation of Springel et al. (2005), which followed the evo- 
lution of 10 billion dark matter particles in the AC DM 
model, to study the properties of more than 10^ haloes 
that were formed. Their study showed that the log- 
normal distribution cannot approximate well the tails 
of the resulting distribution so they proposed other al- 
ternatives. An alternative of similar form, but rather 
more general, is used in the present paper. This is of 
the form: 



0(A') 



A 



exp 



(15) 



This distribution peaks at X'peak = •^o(^) ^^"^ 



AX 



pectation value for A' is given by < A' > = 

where A, a,b,c,A'o are positive parameters and F is 
the complete gamma function. In Fig. 3, we present two 



characteristic snapshots for the distributions of A', for 
the cases CI and C5 that correspond to the smallest and 
the largest haloes studied. Squares correspond to the 
results predicted by the method described in this paper, 
dashed lines correspond to the best least-squared fits of 
the results by a log-normal distribution while solid lines 
are the best least-squared fits by a distribution of the 
form of eq.l5. It is clear that distributions given by 
both cq.l4 and eq.l5 are good fits of our results for the 
low mass case CI. For case C5, the log- normal distri- 
bution cannot fit well the tails of the solid line. In both 
cases, the results are described better by eq.l5. It is 
also clear that the distribution depends on the mass of 
the halo. Although the peaks arc located at about the 
same position in both snapshots, the value of the peak 
is significantly larger in the case CI. We also note that 
a\i is a decreasing function of toq. It varies from 0.71 
for case CI to 0.67 for case C5. It has also been noticed 
by other studies (e.g. Vitvitska et al. 2002) that < A' > 
does not depend on the mass of the halo. Fig. 3 reveals 
that this is not true in our results. This can be seen 
more clearly in Fig. 4, where < A' > versus the present- 
day mass of the halo, niQ, is plotted. Although for a 
factor of masses 10000 -from 0.01 to 100- < A' > varies 
by a factor only of about 1.41 - from 0.0343 to 0.0484- it 
is obvious that < A' > is an increasing function of mass. 
Checking the ability of formula (15) to fit the results of 
Monte-Carlo predictions we found that for Aq = 0.049 
and c = 1 we can predict very good fits to all cases, 
where the values of the rest of the parameters are given 
in Table 1 
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Fig. 5 The role of recent major mergers in the distribution of A' for the cases C1,C3 and C5. The solid line in every 
snapshot shows the distribution for all the haloes in every case. Dashed lines show the distribution of the haloes that had 
at least one recent major merger while dotted lines correspond to the haloes that had no recent major merger. 



A and (3 seem to be decreasing functions of mass. As 
regards a, this is clearly an increasing function of the 
mass of the halo. 

According to the results of N-body simulations it is 
likely that spin parameter is a decreasing function of 
mass. This is supported, for example, by the results 
of Bett et al. 2007 and Maccio et al. 2007. Addi- 
tionally, Bett et al. showed that the values of spin 
parameter and its behavior as a function of mass de- 
pends crucially on the halo-finding algorithm. This 
conclusion was derived by studying three different halo- 
finding algorithms, the traditional " friends-of-friends" 
(FOF) algorithm of Davis et al. (1985), the "spherical 
overdensity" SO) of Lacey & Cole (1994) and a new 
halo definition that they introduced, the TREE haloes. 
Since the results are so sensitive to the halo-finding al- 
gorithm there is a problem regarding the comparison of 
the results of N-body simulations that have used dif- 
ferent halo-finding algorithms. Especially, this problem 
is more serious for the results that are relative to the 
growth of angular momentum. In any case, it seems 
that N-body simulations favor a spin parameter that is 
a decreasing function of the virial mass. Instead our re- 
sults as those of Mailer et al. (2002) favor the trend that 
the spin parameter is an increasing function of the virial 
mass of the halo. Taking into account that our merger- 
tree algorithm has been constructed independently of 
that of the above authors, the disagreement with the 
trend seen in N-body simulations is not likely to arise 
from inconsistencies in the merger-tree construction al- 



Table 1 



Parameters 


A 


a 


/? 


Case 1 


272.206 


1.307 


3.302 


Case 2 


222.880 


1.292 


3.032 


Case 3 


183.043 


1.326 


2.76 


Case 4 


178.417 


1.473 


2.67 


Case 5 


187.216 


1.719 


2.66 



gorithm. Additionally the ellipsoidal collapse model we 
use, that is an improvement of the spherical collapse 
model used by the above authors, is also incapable to 
resolve the problem. We also have to note here that a 
similar trend as that seen in our results is also found, 
for A, in the results of Einsenstein & Loeb (1995) where 
the collapse of homogeneous ellipsoids in the tidal field 
of their environment is studied. 

In order to study the role of the distribution of veloc- 
ities we also used a different distribution according to 
the predictions of N-body simulations of Vitvitska et 
al 2002 that is defined as follows: Let that haloes with 
masses M and m and radii R and r merge. For con- 
venience capital letters correspond to the larger halo. 
After the choice v^ei: as described in section 2., we 
find v,„,M = [GM/Rf/^ and v,^r,m = [Gmlrfl^. 
Then the tangential velocity is picked by a Gaus- 
sian with mean value vt^mean ~ (0.9 — 0.5 ^°"''^ )vcir,M 

for Vcir,m/Vcir,M > 0.4 and Vt^mean = 0-7Vcir,AI for 

/vcir,M < 0.4. The above Gaussian has a = 
vt,meanli- This schcmc is consistent to the results of 
Colin et al. (1999) that major mergers are significantly 
more radial than minor ones and bring in less specific 
angular momentum. The differences between the re- 
sults of this new distribution are those derived by the 
distribution described in section 2 are negligible and 
the spin parameter is again an increasing function of 
the virial mass. However the question remains open. 
A number of authors has stressed the role of recent 
major merger in the final value of the spin parameter. 
More specifically, they have shown that haloes which 
have suffered large recent major mergers appear to have 
larger values of spin parameter. We define as recent 
mergers those occurred at redshifts z < 3. We con- 
sider the merge between a large halo M and a small 
m as major if m/M > 0.5. Obviously a recent major 
merger satisfies both the above conditions. Figs. 5 and 
6 show the dependence of the < A' > distribution on 
the number of recent major mergers for three of our 
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cases namely CI, C3 and C5. In Fig. 6 solid lines show 
the distribution of < A' > for all haloes of the respec- 
tive case, while dashed lines are the distributions over 
those haloes that had at least one recent major merger. 
Dotted lines are the distributions of haloes that had no 
recent major mergers. It is clear, particularly at low 
masses as in case CI, that dashed lines represent dis- 
tributions that are shifted to the right relative to that 
represented by solid lines. Thus, haloes that had at 
least one recent major merger have larger spin param- 
eters. On the other hand, haloes that had no recent 
major merger -represented by dotted lines- show a nar- 
row distribution shifted to the left, relative to the solid 
line, that has an obviously smaller mean value. 

In Fig. 6 we present < A' > for various groups of 
haloes from the cases CI and C5 versus the fraction of 
the total number of haloes of the ease that these groups 
represent. Squares correspond to the case CI, stars to 
the case C3 and triangles to the case C5. From the 
left to the right symbols mark the mean value of spin 
parameter versus the fraction for the following groups: 
i) all haloes of the case, ii) haloes that had at least one 
recent major merger, iii) haloes that had at least two 
recent major mergers and iv) haloes that had at least 
three recent major mergers. The square inside the circle 
marks < A' > for those haloes of the case CI that had 
no recent major mergers. The same is indicated by the 
star inside the circle and the triangle inside the circle 
for the cases C3 and C5, respectively. This fig. shows 
that < A' > is an increasing function of the number of 
recent major mergers. For example, the group, from 
case CI, of haloes that had at least one recent major 
merger has < A' >= 0.0461 and represents 35.6% of the 
total haloes of the case while the group of haloes that 
had at least three major merger has < A' 0.0545 
and represents only 4% of the total number of haloes 
of the case. For the group of haloes that had no re- 
cent major merger < A' >= 0.0182, a significantly low 
value, while the haloes of this group represent 24.3% 
of the total number of haloes of the case. It is noticed 
that as the haloes belonging to that last group had an 
unperturbed recent history, they have time to evolve 
their gas smoothly to a rotationally supporting disk. 
Its fraction is an decreasing function of the final halo 
mass. For case C5 only 0.63% belong to that group. 
The fraction of haloes with no recent major merger is 
definitely a decreasing function of the present day mass 
of the halo. It is quite reasonable, in the hierarchi- 
cal clustering scenario studied here, that recent major 
mergers become rare effects for small haloes. It would 
be interesting to see if, at fixed number of recent major 
mergers, the distribution of A' is independent of mass. 
For this reason, in Fig. 7 we plotted the mean value of 



A' for haloes that have suffered the same number of 
recent major mergers(0,l,2,3 and 4) for all cases. Dif- 
ferent symbols correspond to haloes of different masses. 
Squares correspond to CI, stars to C2, diamonds to C3, 
circles to C4 and deltas to C5. It is clear from this fig. 
that for haloes that have suffered the same number of 
recent major mergers, the heavier one has the larger 
A'. Summarizing the results it is yielded that: a) In all 
cases haloes that had at least one recent major merger 
have < A' > larger than those haloes that had none, 
b) The fraction of haloes that had at least one recent 
major merger is larger than the fraction of haloes that 
had no recent major mergers at all. Thus, we can draw 
the following conclusions: 

1) Disk galaxies are found preferentially in small haloes. 

2) Between haloes of the same mass those that host el- 
liptical galaxies rotate faster than haloes of spiral galax- 
ies (see Vitvitska et al. 2002). 

3) The number of haloes that had no recent major 
merges is significantly smaller that the number of haloes 
that had at least one major merger. If we took into ac- 
count that major mergers destroy galactic disks and 
produce spheroidal stellar systems (see, e.g., Barnes 
1999) small mergers probably do not (see, e.g.. Walker, 
Mihos $ Hernquist 1996) it is natural to expect that 
haloes which have suffered at least one major merger 
could not a host a spiral galaxy. However spheroidal 
stellar systems should be more common objects than 
spiral galaxies. 

Fig. 8 depicts the dependence of the distribution of spin 
parameter on the redshift z. Solid lines correspond to 
present day z = while dashed and dotted lines to 
2=1 and to z = 3 respectively. From the left to the 
right, snapshots correspond to the cases CI, C3 and 
C5 respectively. Differences between z = and z ~ \ 
are small. Curves at z = 1 appear with smaller peaks 
and a slight shift, relative to curve for z — 0, to the 
right for small haloes and to the left for larger haloes. 
Differences between z = and z = 3 are more obvi- 
ous. Dashed curves are, for all cases, shifted to the 
right relative to solid curves. The form of the distri- 
butions remains practically unchanged. A shift of the 
distribution to the right shows that the mean value of 
the spin parameter becomes larger while a shift to the 
left shows that it becomes smaller. Thus, the curves 
indicate that the value of the spin parameter decreases 
from z = 3 to z = for all cases. This is verified by 
the straightforward evaluation of < A' > that appears 
to be systematically smaller at z = than its value at 
z = 3. 
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Fig. 8 The dependence of the spin distribution on the redshift z. Solid lines correspond to present day, z = while dashed 
and dotted lines correspond to 2 = 1 and to z = 3. From the left to the right, snapshots correspond to the cases CI, C3 and 
C5, respectively. It is shown that the distribution of A' is, to a good approximation, unchanged between z = G and z = 1. 
Differences between z = Q and z = 2> are more obvious but without a significant change of the shape of the distribution. 



4 Discussion 

This study describes a picture for the growth of the an- 
gular momentum of dark matter haloes in terms of a 
hierarchical clustering scenario. The results presented 
above are, in general, in good agreement with the re- 
sults that were already known in the literature. Com- 
paring our results with those of large N-body simula- 
tions, we have found satisfactory agreement in the fol- 
lowing points: 

1) The values of spin parameter are in the range of 
0.0343 to 0.0484 for haloes with present day masses in 
the range of lO^h-^Mg to IQ^^hT^MQ. 

2) A log-normal distribution approximates satisfacto- 
rily the distributions of the values of the spin parameter 
but it fails to describe accurately the tails of the result- 
ing distributions. A new, more satisfactory formula, is 
presented. 

3) The role of recent major mergers is very important. 
The distribution of the spin parameter is appreciably 
affected by the number of recent major merger. The 
present day value of the spin parameter of a halo is an 
increasing function of the number of the recent major 
mergers. 

4) The distributions of the spin parameter do not de- 
pend significantly on the redshift. 

5) The value of the spin parameter is a function of the 
present day mass of the halo. The form of this function 
depends, in N-body simulations, on the halo-finding al- 
gorithm but in general seems that spin parameter is 
a decreasing function of mass. Instead, in our results, 
< A' > is an increasing function of mass, approximately 
very closely a power-law form. 

Our results give rise to some questions, as for example: 
Why semi-analytical methods are not able to predict 
the correct relation between the spin parameter and the 



virial mass of the halo? Does this disagreement reflects 
the null role of tidal fields in the orbital-merger picture 
or it arises from other problems associated with the na- 
ture of merger-trees? Are merger-trees able to give the 
correct relation for better description of the velocity 
field during the merge? Is any way improvements on 
both analytical and numerical methods are required in 
order to help us answering some of the above questions 
and to advance our understanding about the physical 
processes that created the structures we observe. 
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Fig. 6 The role of recent major mergers for the cases C1,C3 
and C5. The horizontal axis shows values of < A' > while 
the vertical one represent the fractions of various groups 
of haloes. Squares correspond to the case CI, stars to the 
case C3 and triangles to the case C5. From top to bottom 
symbols mark the mean value of spin parameter versus the 
fraction for the following groups: i) all haloes of the case, ii) 
haloes that had at least one recent major merger, iii) haloes 
that had at least two recent major mergers and iii), haloes 
that had at least three recent major mergers. The square in 
the circle marks < A' > for those haloes of the case CI that 
had no recent major mergers. The same hold for the star 
inside the circle and the triangle inside the circle but for the 
cases C3 and C5 respectively. See text for more details. 



Fig. 4 The mean value of A' as a function of the present 
day mass mo of the halo. Although the values of A' are 
not very sensitive to the values of the mass, as they vary 
from 0.0343 to 0.0489 for a range of mass 0.01 to 100, our 
results show that < A' > is an increasing function of the 
mass of the halo. The linear dependence shown, log{< A' > 
) ~ —1.392 + 0.0387Zog(?no), in this figures indicates that 
< A' > varies as a power of mo. 




1 2 3 

Number of major mergers 



Fig. 7 This figure presents the distribution of spin pa- 
rameter at fixed number of recent major mergers (0,1,2,3 
and 4). Difi'erent symbols correspond to haloes of different 
masses. Squares correspond to CI, stars to C2, diamonds 
to C3, circles to C4 and deltas to C5. It is clear that even 
at fixed number of recent major mergers the distribution of 
A' is an increasing function of the mass of the halo. 
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